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Abstract 

We estimate the probability that a disoriented chiral condensate 
forms during the spherical expansion of a hot medium described by the 
linear sigma model. 


1 Introduction 

The suggestion [0-0 that disoriented chiral condensates (DCC) could form in 
high-energy hadronic collisions has attracted much attention. More than hun¬ 
dred publications are devoted to the study of various aspects of this hypotheti¬ 
cal phenomenon (see e.g. the reviews [j|, 0). The most plausible mechanism of 
DCC formation has been identified: fast expansion of hot quark-gluon plasma 
results in a rapid suppression of thermal fluctuations (quenching), which in 
turn triggers a dramatic amplification of soft pion modes (see e.g. refs. [000. 
However, nobody made a serious attempt to estimate the probability that this 
happens. It is true that the present models are not realistic enough to be 
trusted at the quantitative level. Nevertheless, even within the existing frame¬ 
work it is legitimate to seek for a crude estimate of the probability in question. 
This is the problem we address in the present paper. 

It is expected that the chiral symmetry is restored in a sufficiently hot 
quark-gluon plasma and that it is spontaneously broken when the plasma is 
cooled. At the same time, freely propagating quarks and gluons get trapped 
into hadrons. It would be too difficult to model both phenomena, symmetry 
breaking and confinement. Hence, in the DCC literature it is customary to 
describe the soft modes of the cooling medium by the linear sigma model, 
evacuating cavalierly the confinement problem. The justification of this ideal¬ 
ization, which we also adopt, is discussed at length in many places. 

1 Laboratoire associe au Centre National de la Recherche Scientifique - URA00063. 
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The fast expansion can be modeled by assuming that the field depends 
on time via the variable r = 


the index /j running from 1 to D. 
In the equations of motion the operator d 2 /dt 2 is then replaced by d 2 /dr 2 + 
( D/r)d/dr . In addition to the acceleration term, there also appears a ’’fric¬ 
tion 1 ' term reflecting the decrease of energy in a covolume Since the ’’fric¬ 
tion force” is proportional to D, the larger is D the more efficient is the 
quenching 

In this paper we focus on the spherical geometry, the most favourable 
one for DCC formation. We imagine a droplet of DCC undergoing a radial 
expansion in its rest frame. As remarked in |L0| . 


where the same geometry is 
studied, this does not necessarily mean that the full collision process has the 
same symmetry. 

The time evolution of the chiral phase transition during a spherical ex¬ 
pansion of a medium described by the linear sigma model was studied in ref. 
0. We adopt most of their formalism here and, in particular, the two basic 
ingredients: the assumption that the system is initially in local thermal equi¬ 
librium and the self-consistent approximation of the mean-field type. We shall 
argue that essentially no other assumptions are needed to determine the prob¬ 
ability that soft modes are amplified by a given amount. This is, however, not 
quite enough to fix the probability of observing a DCC in a collision process. 
Information on the environement of the DCC bubble and on the experimental 
set-up is needed for that. We shall come back to this question towards the 
end of this paper. 

In the next section we review briefly the formalism. In sect. 3 we explain 
our sampling strategy for initial conditions. The results and discussion are in 
sect. 4. 


2 The formalism 

Let us briefly summarize the formalism of ref. |8j . The authors start with the 
linear sigma model for an TV-component scalar field (eventually N is set to 4). 
Using the 1/iV approximation they derive the following dynamical equations 
of the Hartree type: 

(d 2 + x)4>(x) = Hn a , ( 1 ) 

(< 9 2 + x)4>(x) = 0 , ( 2 ) 

where x satisfies the gap equation 

X = A 0 (0 2 + (0 2 >-u o 2 ), (3) 

<fi is the classical mean-field and <fi is the quantum fluctuation. The parameters 
Aq and r; 0 are the bare couplings of the sigma model. An explicit symmetry 
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breaking term, proportional to H , gives a finite mass to the Goldstone boson. 
Of course, (4> 2 ) diverges and it is necessary to set an ultra-violet cut-off A. 
The quadratic divergence is removed by subtracting the mass gap at zero 
temperature, T = 0, obtaining 

T-(x - m l ) = 0 2 - fl + (0 2 ) - (0 2 )t=0 (4) 

'-'O 

The remaining logarithmic divergence is eliminated by introducing the renor¬ 
malized coupling constant: 

K 1 — K 1 + (N/8n 2 ) [ k 2 dk/ ( k 2 + m 2 ) 3/2 (5) 

Jo 

The cancellation of divergences is, at any time, insured by an appropriate 
quantization condition (cf eqs. (|TT| ) below). The theory becomes trivial when 
A —> oo and therefore the cut-off has to be kept finite. The effect of the 
renormalization is to reduce the cut-off dependence of the results. 

Notice, that in this approximation, unlike in the standard Hartree one, 
all components of the fluctuation field have the same mass y/x, whose value 
is independent of the orientation of the order parameter cf). This is not very 
realistic but is acceptable if one wishes to estimate only the amplification of 
the length of the field vector, without paying attention to its orientation in 
the internal space. 

In order to take care of the expansion, the time t and the radial coordinate 
r are replaced by the hyperbolic coordinates 

r = Vt 2 — r 2 , q = tanh -1 (r/t) (6) 


and the theory is quantized on hypersurfaces r = const. It is assumed that 
the mean-field depends on r only. Projecting the fluctuation field on the 
orthonormal eigenmodes y s (r],d,Lp) of the curved Laplacian (cf Jll|; we use 
the shorthand notation s = ( s,l,m ), s is dimensionless ) one gets 


0s (t) = 0 s (r)a s + (-l) m 0 s (r)al s (7) 

0s(t) = ^ s (r)a s + (-l) m 0 a (r)a! s (8) 


where 

[ajs,Oj s? ] 5ij5 gg? (9) 

and the dot represents the differentiation with respect to r. The problem can 
be reduced to the familiar one of a set of parametrically excited harmonic 
oscillators. For that purpose one replaces if; s —> g s /r and one introduces the 
’’time” variable u = log (m n r) to get 

[^2 +u s (u)]g s (u) = 0 ( 10 ) 
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where u s (u) = \Js 2 3 + x(t)t 2 is the dimensionless ”frequency”. The oscillators 
are coupled through the common mass. The physical frequency is uj s (u)/t. 

The quantization is completed by choosing appropriate initial conditions 
for the mode functions. The cancellation of divergences required by the renor¬ 
malization of the mass gap is insured if one adopts the following adiabatic 
condition^ at u = u 0 : 

g s (u 0 ) = gi°\u 0 ), g' s (u 0 ) = g^iuo) 

with 

gfHu) = 1 l<j2u,(u) (12) 

st’W = -K(“)/ 2w .(“) + i w .(“)bl 0) (“) (13) 

the prime representing the differentiation with respect to u. 

Assume that Heisenberg and Schroedinger representations coincide at 
t = tq. One can then regard a§ as the Schroedinger representation operator 
creating a particle with frequency u s {uq)/tq. These are the particles appro¬ 
priate for the description of the initial state. One can introduce analogously 
the particles with frequencies uj s {uf)/rf appropriate to the final state and 
counted at r = 77 . Let denote the corresponding creation operator, in the 
Schroedinger representation. The Heisenberg representation held operators, 
given by (7)-(||), can then also be written 

0s(t) = ^ 0) (r/)b s (r) + (-l) m ^ 0) (r / )b! s (r) (14) 

</>s (r) = ^ 1) (r/)b s (T) + (-l) m ^ 1) (r / )bl s (T) (15) 

where ^°\r) = gf\u)/r and ^\t) = [(/^(u) — gf\u)\/r 2 while b s (r) = 
U(t, ro)b s f/^(r, To), U being the unitary operator connecting the two repre¬ 
sentations (when the potential is quadratic U can be explicitly constructed 
, but we do not need this construction here0). Using (7)-(§) and (14)- (|T5 ) 
one easily finds the Bogoliubov transformation connecting the operators as, a§ 
and b s (r), bg(r): 

b s (r) = a(r)a s + /3(r)(-l) m al s (16) 

where 

®s{t) = i[g' s (u)gi°\u f ) - g s (u)g^\u f )} (17) 

Ps(r) = i\g' s (u)gi°\u f ) - g s (u)g^(u f )] (18) 

2 Only second order adiabatic condition is imposed. This may be insufficient in the general 
context 0 but is enough when one only seeks to insure the proper renormalization of the 
mass gap. 

3 The operator U combines free propagation and squeezing. The latter takes care of the 
creation of particle pairs. 



( 11 ) 
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Assuming 


< a iS a js’) = n i“ ) ( r o)Mss’ ’ («*sa jS ’) = 0 (19) 

one finally obtain^] , for one field component, the number of b-particles at 
time r for a given number of a-particles at time To(] : 

nS(r) + i = A s (T)[n^ (r 0 ) + (20) 

where 

A s = 2\ p s \ 2 +1 (21) 

The multiplicity ny( r ) calculated at r = Tf depends weakly on the choice of 
the reference final time Tf , when the latter is large. Assuming local thermal 
equilibrium in the initial state one sets n^(r 0 ) = {exp [uj s (u 0 )/t 0 T] — l} -1 . 


3 Physical picture and the sampling strategy 

As already mentioned in the Introduction, we consider the evolution of a spher¬ 
ical droplet of DCC in its rest frame. We start with a small ball of radius R , 
filled with hot matter in local thermal equilibrium. We assume that the ball 
expands at the speed of light. Due to the time dilation the temperature and 
the value of the mean-field stay approximately constant within a layer near the 
boundary of the ball. Of course, the width of this layer shrinks to zero with 
increasing time (this width equals the distance between the surface r = R and 
the light cone). The equations of motion of the sigma model are supposed to 
describe what happens in the interior of the ball, the cooling observed as one 
moves away from the surface towards the center. We make first the unrealistic 
assumption that the DCC is connected forever to a heat bath kept at constant 
temperature, so that the process never stops. Later on we shall discuss the 
effect of switching the heat bath off. 

The evolution of the droplet depends on the initial conditions set at 
r = To = R. It is easy to convince oneself that in the Hartree approximation 
the only relevant initial conditions are those concerning the mean-field and 
its time derivative (the fluctuations of the initial occupancy numbers appear 

4 The derivation given here differs from that of ref. [|j, where the so-called adiabatic 
basis is unnecessarily used. Strictly speaking, the Bogoliubov transformation given by their 
eqs. (3.19) is not unitary when the frequency is imaginary in some time interval, since 
the exponential factor entering the definition of the adiabatic basis is then no longer a pure 
phase. We understand that in actual calculations the amplification factor was obtained from 
the formula derived for real frequencies. We reproduce their results using eqs. (17)-(18). 

°J. Randrup has emphasized that the amplification of vacuum fluctuations contributes 
significantly to the effect, which is therefore likely to be underestimated using classical 
equations of motion |L4| 
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only in the formula for the initial mass gap and approximately average to zero 
there). How to choose these parameters? In all publications up to now they 
were given values ad hoc. We shall argue, that the probability that </>(r 0 ) and 
0 (t o ) take a given value is determined once one has assumed that the droplet 
is initially in local thermal equilibrium. In order to make the point let us for 
a moment neglect the complications due to quantum mechanics and let us us 
consider a classical field. 

Consider a ball of fixed volume Vo filled with the field in local thermal 
equilibrium. Hence, the field fluctuates as if the ball were part of a system 
in equilibrium, with volume V much larger than Vo- In this large system the 
variances of the spacial averages of the field and of its time derivative would 
be very small, of order 1/V (since long range correlations are absent). The 
corresponding variances for averages calculated by integrating over the volume 
of the ball would be larger by a factor of order V/V 0 (assuming that the radius 
of the ball is at least of the order of the correlation length in the medium). 
An observer living within the ball would identify the spatial averages of the 
field and of its time derivative with 0(ro) and </>( to ), respectively. The point 
is that these random variables are fluctuating in a predictable manner. 

In quantum theory 0(ro) and 4>(tq) can be sampled from the probability 


distribution^ characterized by the mean 

E[H To)] = SioH/xx, (22) 

E[k to)] = 0, (23) 

where Xt is the mass gap (squared) in equilibrium at temperature T, by the 
variances 

Var[0i(r o )] = ^ [ d 3 x d 3 x' ((fi(x)(j)i(x')) (24) 

Vq JVq 

Var [$i(T 0 )] = ~2 [ d 3 x d 3 x' (25) 

Jvq 

and by the covariance 

Cov[0i(r o ),0(r o )] = 772 /. d 3 x d 3 x' (^-[^(x^^x') + <j>i(x')(t>i(x)]) (26) 

Vq JVo 2 


In the Hartree approximation the interacting system is replaced by an ensemble 
of free excitations, with mass y/xfr- It is not difficult to convince oneself that in 
this approximation and in the thermal canonical ensemble with density matrix 
oc exp(— Hq/T), where H 0 is the free Hamiltonian, the probability distribution 

®The probability distribution of the measured values of an observable O is determined 
by the characteristic function ( e lko ). 
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we are interested in is Gaussian. Furthermore the covariance vanishes. Hence 
the parameters given by (22)-(25) determine the distribution, actually the 
analogue of the Wigner function, exactly. 

The variances (24)-(25) can be estimated analytically, when the radius R 
of the ball is much larger than the correlation length A = 1 /^/xt- It is sufficient 
to calculate the variances for the quasi-infinite volume V and to multiply the 
result by V/Vq, since the V/Vq small cells fluctuate independently. One gets 


Var [0i(r o )] 
Var [$i{r 0 )] 


3coth(yxF/2 T) 

3y / Xrcoth( v ^r/2T) 
8 vr R 3 


(27) 

(28) 


The dispersion of calculated exactly is smaller by a factor of 2 (1.5) for 
R = A (R = 2A) and T = 200 to 400 MeV. The dispersion of 0; obtained from 
(^8j) is off by 20% (8% ), respectively. For R < A the disrepancy between the 
analytical formulae and the exact results increases rapidly. Thus, as expected, 
the fluctuations within the ball depend rather weakly on the environement 
provided R > A. 

The formalism of ref. [Q, reviewed in sect. 2, together with the sampling 
method proposed in this section enable one to estimate the likelihood of a 
coherent amplification of the pion held. More precisely, one can calculate the 
probability that the amplification factor A s given by eq. (21) takes a given 
value. In such a calculation the size of the initial ball, viz. R , appears as a 
free parameter. Remember, however, that one has to set r 0 = R and that 
the friction force responsible for the quench is proportional to 1/r. Thus 
the likelihood of DCC formation decreases with increasing R: this parameter 
should be assigned the smallest possible value in order to get the estimate we 
are looking for. 

At this point of the discussion let us remark that the theory we wish 
to use makes only sense for R > A. Indeed, the concept of local thermal 
equilibrium is meaningfull when applied to a cell whose degrees of freedom 
fluctuate more or less independently from what happens in the neighbor cells. 
Also, the validity of the Hartree approximation requires the size of the cell to be 
larger than the Compton wavelength of an excitation. With these arguments 
in mind we focus on the values of R in the range of one to two correlation 
lengths []. 


7 We cannot exclude that a DCC instability develops in a smaller cell, but we have no 
theory to deal with such a scenario. 
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4 Results and discussion 


We show in Fig. 1 the histograms of the probability P(A) that the amplifica¬ 
tion factor of the s = 0 + mode is A 0 = A . It was convenient to assign to the 
parameters of the model the values already used in ref. H[| Proceeding in this 
way we could check our results against theirs. The amplification is calculated 
at r = Tf = lOfm, where the system is in the stationary regime. It is seen 
from the figure that P(A) falls rapidly with A. Clearly, large amplification 
occurs in a small fraction of events only, especially with the choice R = 2A. 
This is qualitatively similar to the analytical result obtained in ref. 0 using a 
1+1 dimensional toy model, except that in the present case P(A) has a power 
falling tail at large A: 

The histogram corresponding to T = 200 MeV and R — A is fairly well 
represented by 


P(A) 


0.805 

(1 + 0.27 A) 3 - 21 


(29) 


The corresponding fit for T = 400 MeV is 


P(A) 


0.213 

(1 + 0.032k!) 6 - 959 


(30) 


It is impossible to find a good fit reproducing the tail of the histogram with 
an exponential function, e.g. with P(A) oc exp (— aA b ). 

It is interesting to inquire what characterizes the initial states leading 
to large amplifications. It appears that one condition is the smallness of the 
absolute strength of the initial classical isovector current (calculated using 
0 (to) and This is ilustrated in Fig. 2, where we show the square of 

the isovector (resp. isoaxialvector) current versus the amplification factor A. 

In order to judge what amplification should be considered as ’’large” one 
has to estimate the multiplicity of produced pions. This can be done with the 
help of the formula (cf |H|) : 


cLyr r 

E dh = I d X '/-Q S ( T ~ r /)/b.P)P'X 


(31) 


where f(x,p) is the invariant phase-space density and u M = x^/r is a unit 4- 
vector orthogonal to the hypersurface r — Tf, where the particles are counted^. 
Eq. ( |3T| ) becomes identical to the eq. (B2) used for the same purpose in ref. 

s Hence, A = 800 MeV, X r = 7.3, /„• = 92.5 MeV. We find the correlation length A = 1.17 
fm (0.68 fm) at T = 200 Mev (400 MeV). 

9 The textbook formula for the intensity of the black-body radiation is obtained replacing 
in d3l| ) the constraint r = t/ by t =const, so that = E, and by setting f(x,p) = 

2(27r)~ 3 [exp (E/T) — l] -1 . The factor 2 is the number of photon polarization states. This 
example helps fixing the normalizations. 







Figure 1: The amplification factor A of the softest mode for (a) T = 200MeV and 
R = A (3.2 x 10 4 MC events) , (b) T = 400MeV and R = A (10 4 MC events) , (c) 

T = 200MeV and R = 2A (2.5 x 10 4 MC events) , (d) T = 400MeV and R = 2A 

(10 4 MC events) . 

||], when one goes over to the hyperbolic coordinates (p;) and integrates out 
the delta function^!. We further set 

f(x,p) = N(27r)- 3 n^ ( 32 ) 

and, as in |8|], we relate s to the 4-momentum by the obvious relation 

= \l{s/Tf) 2 + ml ( 33 ) 

The integrand in ( PTTI ) depends on a single external 4-vector, viz. and 

10 Except for a factor 27r, coming from the integration over the azimuthal angle, which is 
missing in (B2) due to a typo. 
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Figure 2: [resp. ], the square of the initial classical isovector [resp. 

isoaxialvector ] current (in fm -6 ) versus the amplification factor A. We use Monte 
Carlo data corresponding to R = A and T = 200 MeV (3.2 x 10 4 MC events). 


therefore the invariant spectrum is flatj 11 : 


„ dn . 

E dF P = cA 


(34) 


The constant c on the rhs is varying from one Monte Carlo event to 
another, but for A > 30 its average is roughly c = 5GeV -2 for R = A and both 
T = 200 MeV and T = 400 MeV. The flatness of the spectrum is, of course, 
an artifact of the unrealistic assumption that the boost invariant expansion 


continues forever. It is worth mentioning that the rhs of (|34|) does not depend 
on the choice of 7y, provided the input f(x,p ) in ( |3T| ) is a time independent 
function of s. 


ii 


We are puzzled by the falling DCC spectrum shown in fig. 7 of ref. 
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At this point one should distinguish between the intrisic DCC dynamics 
and the extrinsic aspects of DCC formation, those determined by the behavior 
of the environement of the DCC bubble. A discussion of the latter, which 
would require the modeling of the collision proces as a whole, is beyond the 
scope of this paper. Let us only remark that in a real collision process the 
expansion would last a finite time. The resulting spectrum would be cut, the 
value of the cut-off reflecting the behavior of the environement. Unfortunately, 
the predicted total multiplicity depends strongly on this cut-off and cannot 
be estimated in a reliable manner. However, the rhs of (|34[) is presumably 
a reasonable estimate of the invariant momentum space density of soft pion 
radiation. 

A simple example is instructive: Denoting by h the height of the rapidity 
’’plateau”, one can parametrize the one-particle spectrum in the central rapid- 
dn — h ‘ Pt ^ In a central Pb-Pb colision at the 

per unit rapidity, i.e. for all pions 


ity region writing E exp (—gy) 


161 about 200 n 


CERN SPS one observes 
h ~ 600. Thus the invariant momentum space density at very small transverse 
momentum is roughly 1900 GeV~ 2 , where we have used (p 2 ) = O.lGeV 2 . The 
corresponding density fluctuation is expected to be of the order of 75 GeV -2 . 
The rhs of (|34| ) should be significantly larger than the fluctuation, for the DCC 
signal to be detectable^. The signal would be more than three times the fluc¬ 
tuation if A > 45. Setting R = A, the corresponding probability is roughly 
4 x ICE 3 . Of course, this is a conditional probability, as it has been assumed 
that the initial plasma droplet was formed in the collision^. 

Summarizing: The mean-held approximation of ref.|| is combined with 
a specific sampling method designed to generate the initial values taken by 
the mean-held and its time derivative. The method rests on the assumption 
that the medium is initially in a state of local thermal equilibrium. The 
long-wavelength excitations of the medium are modeled with the linear sigma 
model. The probability that the amplification of the soft modes has a given 
magnitude has been estimated. The probability of an observable DCC signal 
appears small. A crude estimate indicates that in a central Pb-Pb collision at 
CERN SPS this probability is at best of the order of ICR 3 . This indication 
should be taken into account by experimenters designing DCC hunt strategies. 

Acknowledgements: We acknowledge useful conversations/correspon¬ 
dence with F. Cooper, B. Jancovici, J. Jurkiewicz, M.A. Lampert, A.H. Muel¬ 
ler, J.Y. Ollitraut, R. Omnes and J. Randrup. 


12 Hopefully, a DCC signal can be distinguished from a large fluctuation due to its specific 
charge structure. 

13 We get almost the same estimate at 200 and 400 MeV. Although at 200 MeV the initial 
fluctuations are smaller, the friction is weaker, since we take a larger initial ball. 
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